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Abstract 


Model Based Parameter Estimation (MBPE) is presented in conjunction with the hybrid 
Finite Element Method (FEM)/Method of Moments (MoM) technique for fast computation of the 
input characteristics of cavity-backed aperture antennas over a frequency range. The hybrid FEM/ 
MoM technique is used to form an integro-partial-differential equation to compute the electric 
field distribution of a cavity-backed aperture antenna. In MBPE, the electric field is expanded in a 
rational function of two polynomials. The coefficients of the rational function are obtained using 
the frequency derivatives of the integro-partial-differential equation formed by the hybrid FEM/ 
MoM technique. Using the rational function approximation, the electric field is obtained over a 
frequency range. Using the electric field at different frequencies, the input characteristics of the 
antenna are obtained over a wide frequency range. Numerical results for an open coaxial line, 
probe-fed coaxial cavity and cavity-backed microstrip patch antennas are presented. Good 
agreement between MBPE and the solutions over individual frequencies is observed. 
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List of Symbols 


V Del operator 

V' Del operator over the source coordinates 

e r Dielectric permittivity of the medium in the cavity 

e rc Dielectric permittivity of the medium in the coaxial feed line 

5 ?0 Kronecker delta defined in equation (22) 

Dielectric permeability of the medium in the cavity 
p p -coordinate of the cylindrical coordinate system 

p Unit normal vector along the p -axis 

0) Angular frequency 

AWE Asymptotic Waveform Evaluation 

MBPE Model Based Parameter Estimation 

A^ q \k 0 ) g th derivative of A(k ) with respect to k; ~^—A(k) , evaluated at k 

dk q 

a n Coefficients of the numerator of the rational function (n = 0,1,2, 3, , L) 

B(k ) Excitation vector 

B (q \k 0 ) q lh derivative of B(k) with respect to k ; ~^—B(k) , evaluated at k 

dk q 

b m Coefficients of the numerator of the rational function (m = 0, 1, 2, 3, , M ) 

C r s Binomial coefficient 

ds Surface integration with respect to observation coordinates 

ds' Surface integration with respect to source coordinates 

E Electric field 

E inp Electric field at the input plane S- 

e(k) Electric field coefficient vector 
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Incident electric field due to coaxial line at the surface S inp 
Reflected electric field into the coaxial line at the surface S jnp 
Magnetic field at the surface S ap 
Magnetic field at the surface S inp 
Frequency 

Wavenumber at any frequency / 

Wavenumber at frequency f Q 
Magnetic current at the surface S ap 
Normal unit vector 
Polynomial of order L 
Polynomial of order M 
Factorial of number q 

Reflection coefficient at the input plane S inp 

Distance between the source point and the observation point 

Radius of inner conductor of the coaxial feed line 

Radius of outer conductor of the coaxial feed line 

Vector testing function 

Vector testing function at the surface S ap 

Normalized input admittance of the antenna 

Unit normal along Z-axis 
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1. Introduction 


Cavity-backed aperture antennas are very popular in aerospace applications due to their 
conformal nature. Hybrid techniques have become attractive for numerical analysis of these type 
of problems due to their ability to handle arbitrary shape of the cavity and complex materials that 
may be required for the antenna design. The combined Finite Element Method (FEM) and 
Method of Moments (MoM) technique in particular has been used to analyze various cavity- 
backed aperture antennas[l,2]. In the combined FEM/MoM technique, FEM is used in the cavity 
volume to compute the electric field, whereas MoM is used to compute the magnetic current at the 
aperture. Using Galerkin’s technique, an integro-partial-differential equation is formed. The 
cavity is divided into tetrahedral elements and the aperture is discretized by triangles. 
Simultaneous equations are generated over the subdomains and are added to form a global matrix 
equation. This results in a partly sparse and partly dense symmetric complex matrix, which can be 
solved either by a direct solver or by an iterative solver. The electric field hence obtained is used 
to compute the radiation characteristics and input characteristics of the antenna. 

In most practical applications, input characteristics such as input impedance or input 
admittance are of interest over a frequency range. To obtain the frequency response of the 
antenna, one has to repeat the above calculations at every incremental frequency over the 
frequency band of interest. If the antenna is highly frequency dependent, one needs to do the 
calculations at fine increments of frequency to get an accurate representation of the frequency 
response. For electrically large cavities with large apertures, this can be computationally intensive 
and in some cases computationally prohibitive. To alleviate the above problems, recently 
Asymptotic Waveform Evaluation (AWE) was applied to frequency domain electromagnetics to 
obtain the frequency response[3-6]. In AWE, the unknown electric field was expanded in Taylor 
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series around a frequency. The coefficients of Taylor series were obtained using the frequency 
derivatives of the integro-partial-differential equation resulting from the combined FEM/MoM 
technique [6]. 

In this report, a similar but more flexible method called Model Based Parameter 
Estimation (MBPE)[7,8] is applied for predicting the input characteristics of cavity-backed 
aperture antennas over a wide band of frequencies using the combined FEM/MoM technique. In 
MBPE technique, the electric field is expanded as a rational function. The coefficients of the 
rational function are obtained using the frequency data and the frequency derivative data. Once 
the coefficients of the rational function are obtained, the electric field in the cavity can be obtained 
at any frequency within the frequency range. Using the electric field, the input characteristics such 
as the input impedance or admittance can be calculated. If the frequency derivative information is 
known for more than one frequency, a rational function matching all the samples can be obtained 
resulting in a wider frequency response. 

The rest of the report is organized as follows. In section 2, MBPE implementation for the 
combined FEM/MoM technique is described. Numerical results for an open coaxial line, a coaxial 
cavity fed by a coaxial line and cavity-backed microstrip patch antennas are presented in section 
3. The numerical data are compared with the exact solution over the bandwidth. CPU timings are 
given for each example. Concluding remarks on the advantages and disadvantages of MBPE are 
given in section 4. 
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2. MBPE Implementation for the combined FEM/MoM 

Technique 

The geometry of the problem to be analyzed is shown in Figure L For linear, isotropic, 
and source free region, the electric field satisfies the vector wave equation: 

Vx^|- VxE k 2 e r E = 0 (1) 

where ]i r , e r are the relative permeability and relative permittivity of the medium in the cavity. 
The time variation exp(j(or) is assumed and suppressed throughout this report. Applying the 
Galerkin’s technique, equation (1) can be written in “weak form” as [1] 

JJJ (VxT) • frf VxE W* 2 e JJJt . E dv - Mi J J (T x n) . H ap ds 

y ^ 


= J T • (n x H inp )ds (2) 

^inp 

where T is the vector testing function, S ap is the aperture surface, and S inp is the input surface 
(see fig. 1). H ap is the magnetic field at the aperture and H f . is the magnetic field at the input 
surface. 

In accordance with the equivalence principle [9], the fields inside the cavity can be 
decoupled to the fields outside the cavity by closing the aperture with a Perfect Electric Conductor 
(PEC) and introducing the equivalent magnetic current. 

M = Exz (3) 

over the extent of the aperture. Making use of the image theory, the integrals over S ap in equation 
(2) can be written as 
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J J (T x n) • H ap ds 


s M s„„ 


Sap S ap 


( 4 ) 


where T 5 = Txn and R is the distance between source point and the observation point. V' indi- 
cates del operation over the source coordinates and ds' indicates the surface integration over the 
source region. 

Though the analysis presented in this report is not limited to any specific input feed 
structure, we restrict the presentation of the formulation to the coaxial line as the input feed 
structure. The cross section of the coaxial line is shown in Figure 2. Assuming that the incident 
electric field is the transverse electromagnetic (TEM) mode and the reflected field also consists of 
TEM mode only, the electric field at the input plane S inp is given by 

E inp = ^inc eX P (~J k KcZ) + <V eX P U k J^Z) (5) 


where 


e 

tnc 


= P- 


2nln 


r _2\ 

'V 


( 6 ) 


and 


®re/ ^o^inc 


R is the reflection coefficient and is given by 


R 0 = exp j J E»f£\fe-exp(-2/*^z,) 

M) - 


( 7 ) 


( 8 ) 
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r 2 is the outer radius and r, is the inner radius of the coaxial line. e rc is the relative permittivity 
of the coaxial line. 

Using equation (5) to calculate H inp , the surface integral over S inp in equation (2) can be 
written as 

.wj | T • (n x H inp )ds 

Sinp 



x -G>fr*GW 


+ 


2jkji; c exp(-jkji; c z l ) 




( 9 ) 


Substituting equation (4) and (8) in equation (2), the system equations for the combined FEM/ 
MoM technique can be written as 
JJJ — (VxT) • (VxE)Jv’-* 2 £ r JJjT • E dv 

V V 


-fj r V ■ (J I ♦ rJ I < v ‘ T 4l J <*' * 


+ 





2jkJi; c cxp(-jkJ^ c z x ) 




( 10 ) 
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The volume of the cavity is subdivided into small volume tetrahedral elements. The 
electric field is expressed in terms of the edge vector basis functions [10], which enforce the 
divergenceless condition of the electric field explicitly. The vector testing function is also 
expressed in terms of the edge vector basis functions following the Galerkin’s method. The 
discretization of the cavity volume into tetrahedral elements automatically results in discretization 
of the surfaces S a p and S in ^ into triangular elements. The volume and surface integrals in 
equation (10) are carried out over each element to form element matrices and the element 
matrices are assembled to form global matrices. Equation (10) can be written in matrix form as 

A(k)e(k) = B(k) (11) 


A (it) is a partly sparse, partly dense complex symmetric matrix, B(k) is the excitation vector, and 
e(k) is the unknown electric field coefficient vector. A(k) is evaluated as a sum of three matrices. 

A(k) = A^(k) + A 2 (k) + A 3 (&) + A 4 (k ) (12) 

where 

A,(k) = (VxT) • (VxE)dv-k, 2 £ r IJ1 T • E dv (13) 

J V V 




s ap ' S ap 


A,(k) = J <VT.)|j J (V-»M) eX Py* R W j& 


\ 1 E -(p>) 

2 7t ln(r ~/ S ‘" P 1 S " 


(14) 

(15) 

(16) 
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The matrix equation (11) is solved at any specific frequency, f u (with wavenumber k a ) 
either by a direct method or by an iterative method. The solution of the equation (11) gives the 
unknown electric field coefficients which are used to obtain the electric field distribution. Once 


the electric field distribution is known, the input reflection coefficient can be calculated using 
equation (8). The input plane is placed at z } = 0 , and the reflection coefficient is calculated as 



The normalized input admittance at S t is given by 


Y =! zL 

,n 1 +T 


The input admittance given in equation (18) is calculated at one frequency. If one needs 
the input admittance over a frequency range, this calculation is to be repeated at different 
frequency values. Instead MBPE [7,8] can be applied for rapid calculation of input admittance/ 
impedance over a frequency range. MBPE technique involves expanding the unknown coefficient 
vector as a rational polynomial. The coefficients of the rational polynomial are obtained by 
matching the function and its frequency derivatives of the function at one or more frequency 
points. 

The solution of equation (1 1) at any frequency f 0 gives the unknown electric field coefficient 
column vector e(k Q ), where k Q is the free space wavenumber at f Q . Instead e(k) can be written as a 
rational function. 
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where 


e(k) = 


P L (k) 

QmW 


( 20 ) 


2 3 L 

P L (k ) = a 0 + a^ + a 2 k + a 3 k + + a L k 

Q M (k) = b a + b l k + b 2 k 2 + b 3 k 3 + + b M k M 

b a is set to 1 as the rational function can be divided by an arbitrary constant. The coefficients of 
the rational function are obtained by matching the frequency-derivatives of e(k). If equation (14) 
is differentiated t times with respect to k, the resulting equations can be written as [8] 

e Q\! — P L 


e Qm + e Qht - p L 

e"Q M + 2e'Q M ' + eQ M " = P L " 

e'"Q M + 3e"Q M ' + 3e'Q M " + eQ M "' = P L "' 


(0 


Qm 


+ te 


O-D^(l) 


Q 


M 


+ c 


t,t-m 


e Q + 


+ e Q ( M 


= P 


(0 


where C = — is the binomial coefficient. The system of (t+1) equations provides the 

r ' s s!(r-s)! 

information from which the rational function coefficients can be found if t > L + M + 1 . If the 
frequency derivatives are available at only one frequency f 0 , the variable in the rational function 

can be replaced with ( k-k D ) i.e., 


e(k) = 


Q M (k~k 0 ) 


( 21 ) 
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and the derivatives can be evaluated at k = k 0 . The coefficients of the rational function can be 
obtained from the following equations: 


1 

0 

0 

0 

0 


a o = e ( k o ) 


0 0 

~ e i ~e 0 0 

~e 2 -e x 0 

-e-x -e 0 0 


~ e L + M-l e L + M - 2 • ~ e l 


~ 

p -1 


p 


a l 


e l 


a 2 


e 2 


a L 


e L 


v 


e L + M 


( 22 ) 


(23) 


e (m) 

where e m = — - . For example for a rational function with L=5 and M=4, the matrix equation 
can be written as 


1 0 0 0 0 -e 0 0 0 0 


a, 



0 1 0 0 0 -e x -e 0 0 0 


a 2 



0 0 1 0 0 -e 2 —e x -e 0 0 


a 3 



0 0 0 1 0 -e 3 -e 2 -e x -e 0 


a 4 


e 4 

0 0 0 0 1 -e 4 -e 3 -e 2 -e x 


a 5 

= 

e 5 

0 0 0 0 0 - e 5 -e 4 -e 3 -e 2 


h 


e 6 

0 0 0 0 0 -e 6 -e 5 -e 4 e 3 


b 2 


e l 

0 0 0 0 0 -e 7 -e 6 -e 5 -e 4 


b 7> 



0 0 0 0 0 —eg -e 7 —e 6 -e 5 


h 


e 9 


(24) 


This approach is same as the Pade approximation given in [3], This method has been 
successfully applied to electromagnetic scattering from cavity-backed apertures using a hybrid 
finite element and method of moments technique [1 1], 


13 



If the frequency derivatives are known at more than one frequency, then the expansion about 
k=k 0 cannot be used and the system matrix to solve the rational function coefficients takes a 
general form [8]. For the sake of simplicity, only a two frequency model is presented here. 
Assume that at two frequencies ,fj (with free space wavenumber k,) and f 2 (with free space 
wavenumber k 2 ), four derivatives are evaluated at each frequency. Hence 10 samples of data are 
available (two frequency samples and a total of eight frequency derivative samples) to form a 
rational function with L-5 and M—4 


e(k) = 

Equation (25) can be written as 


, .2 .3 .4 ,5 

a 0 + ci\k + + a^k + + a^k 

1 + b l k + b 2 k 2 + fe 3 * 3 + b 4 k 4 


(25) 


(1 +b x k + b 2 k 2 + b 3 k 3 + b A k 4 )e(k) = a 0 + a,* + a 2 k 2 + a 3 k 3 + a 4 k 4 + a 5 k 5 (26) 

Differentiating equation (26) four times at each frequency, the matrix equation for the solution of 
the coefficients of the rational function (equation (25)) can be written as 


1 

*1 

k 2 

K \ 

*? 

*1 

A 

J°h 

e \ K i 

e (0) k 2 
e \ K i 

e (0) * 3 
e \ K i 

e (0) 

e \ 

A 

0 

1 

2 k x 

3 k ] 

4k i 

5 k 4 

A / 27 

^28 

m 29 

m 2 

10 

0 

0 

2 

6*, 

12*, 

20* 3 

m 37 

m 38 

m 39 

m 3 

10 

0 

0 

0 

6 

24 * , 

60 *^ 

47 

^48 

m 49 

m 4 

10 

0 

0 

0 

0 

24 

120*, 

m 51 

M 58 

m 59 

m 5 

10 

1 

k 2 

k\ 

A 

k 2 

k 5 

k 2 

e 2 )k 2 

efk \ 

eft 

e (0] 

e 2 

A 

0 

1 

lk 2 

3k 2 2 

4k\ 

5*2 

m 71 

m 78 

m 19 

m 7 

10 

0 

0 

2 

6k 2 

12*2 

20*2 

m 87 

^88 

m 89 

m 8 

10 

0 

0 

0 

6 

24* 2 

60*2 

m 97 

a/ 98 

m 99 

Mg 

10 

0 

0 

0 

0 

24 

120*2 

i 

WlO 8 

M l0 9 

"10 

10 



(0) 




e \ 




e {1) 

e \ 


a \ 


( 2 ) 




e \ 


a 2 


( 3 ) 


a 3 


e \ 




( 4 ) 


a 4 



e \ 


a 5 


e {0) 

e 2 


b \ 


(1) 




e 2 


b 2 


J2> 


b 3 


e 2 




( 3 ) 


b 4 


e 2 


( 4 ) 

J 

f 2 _ 


(27) 
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,m m 

where *• = and 

dk dk 

M 27 = -(e\ l) k l+ e\ 0) ) 

M 28 = -{e\ X) k\ + 2e ( ^k x ) 

A / 29 = + 

^2 10 = -(e\ X) k\ + Aefk\) 

M 31 = -(e\ 2) k l+ 2e\ ]) ) 

M 3S = -(e^k 2 l+ 4e\\ + 2e^) 

M 3g = -(e5 2) ^ + 6ei 1) ikj + 6e ( j 0) ^ 1 ) 

M 3 jq — — (^] + 8e, k ] + 12e, /rj) 

M 47 = -(^ 1 + 3 ^ 2) ) 

M 48 = -{ e {2) k] + 6e {2 ) k x + 6e { ^) 

M 49 = + 9e (2) k] + Mtfj 0 *, + 6e ( j 0) ) 

M 4 io = + I2e (2) k\ + 36e\ ]) k] + 24e\ 0) k l ) 

M 51 = ~(e\ 4) k i+ 4e\ 3) ) 

M 58 = -{e {4) k\ + %e {2) k x + 12e ( , 2) ) 
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M 59 = -{e (4) k\ + \2e { ?\] + 36e ( ?k x + 24e\ 0) ) 

M 5 10 = -{e (4) k\ + 16 e ( x ] k] + lle^k] + 964% + 24e ( , 0) ) 

Mu = -(4% + 4 0) ) 

m 78 = -( 4 1 } 4 + 2 4 %) 

m 78 = -( 4 ,) 4 + 34 0) 4 ) 

My 10 = —(^2 ^2 + 4e 2 ^ 2 ) 

M g y = -{efk 2 + 24°) 

m 88 = -(4 2) 4 + 4 4% + 2 4 0) ) 

m 89 = -( 4 2) 4 + 64 ,) 4 + 64 %) 

m & io = -(4 2) *2 +8 4 1) 4 + i2 4 0) 4) 

M 97 = -( 4 % + 3 e^ 2) ) 

M 98 = —(4 ^2 ^4 ^2 4 " 6^2 

m 99 = -( 4 3 ) 4 + 9 4 2 ) 4 + 18 4 1 ) ^ 2 + 6 4 0) ) 

M 9 10 = -(4 3) 4 + 124% + 364% + 244%) 

M 10 7 = -(4% + 44 3) ) 

M ,0 8 = -(4 ^2 + 8e 2 ^2 + l 2e 2 ) 


16 



M,o 9 = -{e^kl + 12<4 3) *2 + 36^2 2) k 2 + 24e^ 0) ) 

m,o 10 = -(4 4) 4+ i64 3) 4 +72 4 2) 4 +96 4% +24 4 0) ) 

In the above equations, e ' 1 ’ , the / th derivative is obtained using the recursive relationship, 


« (0 = A-'tk)\ 


B 


<0 -I (i-s^c,/W 


9 ) 


(*) 


( 28 ) 


L g = 0 J 

A {q) {k) is the q th derivative with respect to k of A(k) and B U \k) is the i h derivative with 
respect to k of B(k). The Kronecker delta b qo is defined as 


r 1 9 = ° 

*- ■ { o , * 0 (29) 

The frequency derivatives of A(k) and B(k) are evaluated and given in [6]. 

The above procedure can be generalized for multiple frequencies with 
frequency-derivatives evaluated at each frequency to increase the accuracy of the rational 
function. Alternatively, the two-frequency-four-derivative model can be used with multiple 
frequency windows. As the complexity of the matrix equation to solve for multiple-frequency- 
multiple derivative model increases with the number of frequency points and number of 
derivatives taken at each frequency, the two-frequency-four-derivative model is followed in this 
report. 


3. Numerical Results 

To validate the analysis presented in the previous sections, a few numerical examples are 
considered. Calculation of input characteristics over a frequency range are done for an open 
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coaxial line, coaxial cavity, and cavity-backed square and circular microstrip patch antennas. The 
numerical data obtained using MBPE are compared with the results calculated at each frequency 
using the computer code CBS3DR[12], which implements the combined FEM/MoM 
technique[2]. We will refer to the latter method as “exact solution ” Due to the hybrid FEM/MoM 
technique, matrix A(k c ) is partly sparse and partly dense. The Complex Vector Sparse Solver 
(CVSS) [13] is used to LU factor the matrix A(k a ) once, and the moments are obtained by 
backsolving the equation (28) with multiple right-hand sides. All the computations reported 
below are done on a SGI Indigo2 (with IP22 processor) computer. 

(a) Open Coaxial line: 

An open coaxial line radiating into an infinite ground plane (fig. 3a) is considered. A finite 
length of the line is used for FEM discretization. The input plane S inp is placed at z = 0 plane 
and the radiating aperture at z = lcm . The discretization of the coaxial line resulted in 1 1 19 total 
unknowns, and the order of the dense matrix due to MoM is 144. One-frequency MBPE with L=5 
and M=4 at/ 0 =7.5GHz is used to calculate the frequency response of the input admittance. Two- 
frequency MBPE at /y=6GHz and f2=9 GHz with L=5 and M—4 is also used to calculate the 
frequency response. The frequency response over the frequency range lGHz-13GHz is plotted in 
Figure 3(b) along with the exact solution calculated at 23 discrete frequency points over this 
frequency range. Both one-frequency and two-frequency MBPE frequency responses are 
calculated at 0.1GHz increments. One-frequency MBPE took 97 secs to generate moments, 
whereas two-frequency MBPE took a total of 166 secs to generate the moments at both 
frequencies. The exact solution took 990 secs to calculate input admittance at 23 frequency values 
from 1GHz to 13GHz. It can be seen that one-frequency MBPE agrees well with the exact 
solution over the frequency range 4.5GHz to 13GHz, whereas two-frequency MBPE agrees well 
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with the exact solution over the frequency range 1GHz to 13GHz. Both one-frequency MB PE and 
two-frequency MBPE are faster than the exact solution over the frequency range. Two-frequency 
MBPE has advantage over the one-frequency MBPE as it requires less computer memory 1 . 

(b) Open Coaxial Cavity: 

An open coaxial cavity fed by a 50Q coaxial line (fig. 4) is considered as a second 
example. The input plane S inp is placed at z = 0 plane and the radiating aperture at 
z = 0.952 cm plane. The cavity volume is discretized using tetrahedral elements, which resulted 
in 4541 total unknowns and the order of the dense matrix due to MoM is 666. The frequency 
response of the return loss (=201og|T| ) is calculated using one-frequency MBPE with L=5 and 
M=4 at / 0 =7.5GHz and is plotted in Figure 5. Two-frequency MBPE with L=5 and M-4 at 
/y=6GHz and/ 2 =10GHz is also used to calculate the frequency response of the return loss along 
with the exact solution calculated at 24 frequency points over the frequency range 1GHz to 
13GHz. It can be seen that the one-frequency MBPE agrees well with the exact solution over the 
frequency range 1GHz to 10GHz, whereas the two-frequency MBPE agrees well with the exact 
solution over the frequency range 1GHz to 13GHz. Both one-frequency and two-frequency 
MBPE frequency responses are calculated with 0.1 GHz increments. One-frequency MBPE took 
1563 secs to generate the moments, whereas two-frequency MBPE took a total of 2834 secs to 
generate the moments at both frequencies. The exact solution took a total of 18,576 secs to 
calculate return loss at 24 frequency points in the frequency range 1GHz to 13GHz. Both one- 
frequency MBPE and two-frequency MBPE are faster than the exact solution for the frequency 
response calculation. 


1. Please see the comment on storage at the end of this section for detailed explanation. 
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(c) Cavity-Backed Square Microstrip Patch Antenna: 

A cavity-backed square microstrip antenna radiating into an infinite ground plane (fig. 6) 

is considered. The input plane S inp is placed at z = 0 plane and the radiating aperture at 
z = 0.16cm. The discretization of the cavity volume resulted in 2,160 total unknowns and the 
order of the dense matrix due to MoM is 544. The frequency response of the input impedance 
(1/T ( „) is calculated using one -frequency MBPE with L= 5 and M=4 at/ 0 =4GHz and also using 
two-frequency MBPE with L= 5 and M— 4 at f/= 3GHz and f2 = 5GWz. The numerical data are 
plotted in Figure 7 along with the exact solution calculated at 23 frequency points over the 
frequency range 1GHz to 7GHz. It can be seen from Figure 7 that one-frequency MBPE agrees 
well with the exact solution over the frequency range 1GHz to 6GHz, whereas the two-frequency 
MBPE agrees well with the exact solution over the frequency range 1GHz to 7GHz. Both one- 
frequency and two-frequency MBPE frequency responses are calculated with 0.01 GHz 
increments. One-frequency MBPE took 1 107 secs of CPU time to generate the moments, whereas 
the two-frequency MBPE took a total of 1 120 secs of CPU time to generate moments at both 
frequencies. The exact solution took a total of 1 1,891 secs of CPU time for computations at 23 
frequency points over the frequency range 1GHz to 7GHz. One-frequency MBPE and two- 
frequency MBPE are faster than the exact solution for the frequency response calculations. 

(d) Cavity-Backed Circular Microstrip Patch Antenna: 

A cavity-backed circular microstrip antenna radiating into an infinite ground plane is 

shown in Figure 8. The input plane S is placed at z = 0 plane and the radiating aperture at 
z = 0.16cm. The discretization of the cavity volume resulted in 6,363 total unknowns and the 
order of the dense matrix due to MoM is 469. The frequency response of the input impedance 
(1 /y . ) is calculated using one-frequency MBPE with L= 5 and M-4 at/ 0 =6GHz and also using 


20 



two-frequency MB PE with L=5 and M=4 at / y =5GHz and / 2 =7GHz. The numerical data are 
plotted in Figure 9 along with the exact solution calculated at 17 frequency points over the 
frequency range 3GHz to 9GHz. It can be seen from Figure 9 that one-frequency MBPE and two- 
frequency MBPE agree well with the exact solution over the frequency range 3GHz to 9GHz. 
Both one-frequency and two-frequency MBPE frequency responses are calculated with 0.01 GHz 
increments. One-frequency MBPE took 913 secs of CPU time to generate the moments, whereas 
the two-frequency MBPE took a total of 1636 secs of CPU time to generate moments at both 
frequencies. The exact solution took a total of 7,650 secs of CPU time for computations at 1 8 
frequency points over the frequency range 3GHz to 9GHz. One-frequency MBPE and two- 
frequency MBPE are faster than the exact solution for the frequency response calculations. 


Comment on Storage: In all the above examples, when solving a matrix equation, one needs to 
store a complex, partly sparse and partly dense matrix A(k 0 ) for exact solution at each 
frequency. In one-frequency MBPE one needs to store the derivative matrices ( A^ q \k Q ) , 
q=l,2,3,...(L+M)), along with the matrix A{k 0 ) . For electrically large problems, this could 
impose a burden on computer resources. This problem can be overcome by storing the derivative 
matrices, A^ q \k 0 ) out-of-core, as the derivative matrices are required only for matrix-vector 
multiplication. In two-frequency MBPE, one needs to store only + ■y- — - 1 derivative matrices 

along with the matrix A(k) at each frequency. Once the moments are calculated at one 
frequency, the memory used for the matrices can be reutilized to generate moments at the second 
frequency, hence reducing the burden on computer memory requirements. In all the numerical 
examples presented with L=5 and M=4, one-frequency MBPE had to store 10 matrices, whereas 
two frequency-MBPE had to store only 5 matrices at each frequency. The memory to store the 
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matrices at one frequency is reutilized to store the matrices at the second frequency. Hence, even 
though the CPU timings for two-frequency MBPE is more than the one-frequency MBPE, if 
computer memory is a constraint, it is advisable to use two-frequency MBPE as an alternative to 
one-frequency MBPE. 


4. Concluding Remarks 

The MBPE technique is applied to the hybrid FEM/MoM technique to obtain the 
frequency response of the input characteristics of cavity-backed aperture antennas. The frequency 
response of input characteristics of an open coaxial line, coaxial cavity, square microstrip patch 
antenna, and a circular patch antenna are computed and compared with the exact solution. From 
the numerical examples presented in this work, MBPE technique is found to be superior in terms 
of CPU time to obtain a frequency response. It may be noted that although calculations are done 
in frequency increments of 0.1 GHz or 0.01 GHz for the examples presented, the frequency 
response at even finer frequency increments can also be calculated with a very nominal cost. In 
one-frequency MBPE, the frequency response is valid over a certain frequency range. In two- 
frequency MBPE, the two frequency values have to be chosen so as to get an accurate frequency 
response between the two frequency values. To get a wide frequency response for any problem, 
either one- or two-frequency MBPE models have to be used with different frequency values to 
cover the complete frequency range. To be accurate over all frequency ranges, a reliable error 
criteria should be developed, which can be used to sample the frequency points to apply MBPE 
model. Development of such a sampling criteria will make MBPE a very effective tool for 
computational electromagnetics. 
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INPUT 


Figure 1 Geometry of a cavity backed aperture in finite ground plane. 
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Figure 3 (a) Open coaxial line in an infinite ground plane. Inner radius r 2 =lcm, Outer 
radius r 2 = 1.57cm, e r =1.0 and L=1.0cm 
(b) Normalized input admittance as a function of frequency. 
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Figure 4 Geometry of a coaxial cavity in an infinite ground plane. Outer radius of the coax- 
ial cavity=l”. Inner radius of the coaxial cavity=0.0181” and L=3/8”. The cavity is 
fed by a 50Q coaxial line. 
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Figure 6 Cavity -backed square microstrip patch antenna in an infinite ground plane fed by a 

50£2 coaxial line. 
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Frequency (GHz) 


Figure 7 Normalized input impedance versus frequency of the cavity-backed square micros- 
trip antenna (figure 6). 
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Figure 8 Cavity-backed circular microstrip patch antenna in an infinite ground plane fed by 
a 50Q coaxial line. 
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Frequency (GHz) 


Normalized input impedance versus frequency of the cavity-backed circular 
microstrip antenna (figure 8). 
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